Classifying Rational Densities Using Two One-Dimensional Cellular Automata 
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Given a (finite) string of zeros and ones, we report a way to determine if tlie number of ones is less 
than, greater than, or equal to a prescribed number by applying two sets of cellular automaton rules 
in succession. Thus, we solve the general density classification problem using cellular automaton. 

PACS numbers: 05.70.Fh, 05.50.+q, 89.80.+h 
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I. INTRODUCTION 



Cellular automaton (CA) is an extremely simple lo- 
cal interaction model of natural systems Although 
the rules of CA are simple and local, complicated spatial 
patterns can be formed. Moreover, from the computer 
science point of view, CA can be regarded as a very spe- 
cial kind of Turing machine with no internal memory. In 
fact, tailor-made CA can be used to simulate certain log- 
ical operations ||^. Therefore, it is instructive to know 
the power and limitations of CA in computation. 

Since CA is essentially a special Turing machine with 
no internal memory, it is natural to ask if CA can be used 
to perform certain task that requires a global counter. An 
example is the so-called density classification problem: 
Consider a one-dimensional (finite) chain of sites with 
periodic boundary condition. Each site is either in state 
zero or state one. The problem is to change the state 
of every site to zero if the number of zeros is more than 
that of ones in the one-dimensional chain. Otherwise, 
every site is set to state one. The density classification 
problem is trivial if one can invoke an internal memory 
to count the number of zeros in the chain. On the other 
hand. Land and Belew showed that density classification 
cannot be done perfectly using one-dimensional CA ^ . 

Later on, Fuks pointed out that the density classifica- 
tion problem can be solved if we apply two CA rules in 
succession For a one-dimensional lattice of N sites, 
he first applies the Wolfram's elementary CA rule 184 
(the so-called traffic rule ||) [{N -2)/2\ times. Then, 
he applies the Wolfram's elementary CA rule 232 (the so- 
called majority rule) [(iV — 1)/2J times. The combined 
automaton solves the density classification problem suc- 
cinctly 0. 

With slight modifications, Fuks is able to classify den- 
sity in the form 1/n for some integer n > 2 [Q. Then, 
Fuks went on to ask if it is possible to use combined CA to 



classify an arbitrary density pc- Here, we report a simple 
solution to his problem using a combination of two CAs 
provided that pc is a rational number. But before stat- 
ing our CA rules, let us re-state the density classification 
problem in the most general context. 

Consider a one-dimensional lattice making up of N 
sites with periodic boundary condition. The state of each 
site can either be zero or one. One is given a critical den- 
sity Pc between zero and one. The initial density of ones 
in the system is denoted by p. (That is, there are Np 
sites in the system with states equal one.) Our goal is to 
find a combination of two CA rules such that after going 
through the CA, the final states of the sites on the lattice 
becomes all zeros (and all ones) if p < pc (and p > pc, 
respectively). And in the case of p = pc, the final state 
consists of exactly Npc sites with state equals one.^ 

Since the number of sites N is finite, it suffices for us 
to consider the case when pc is a rational number written 
in the form p/q for some relatively prime integers p and 
q. Furthermore, we can also assume that p and p ^ q 
since classifying pc = and 1 are trivial using CA. Now, 
we report the two CA rules we use to solve the density 
classification problem for rational critical density. 



II. MODIFIED TRAFFIC RULE 

Our generalized traffic rule goes as follows: First, we 
regard a site with state one as being occupied by a car. 
Otherwise, that site is empty. Then a car can move in 
the next timestep to the right by one site if and only 
if (a) its immediate right-neighboring-site is unoccupied; 
and (b) the q — 1 consecutive right-neighboring-sites is 
occupied by at most p—l cars. Otherwise, the car stays 
in its original position. For example, if pc = 2/5, then 
the first cars in the sequences 10000 and 10100 can move 
one step forward in the next timestep, while the first cars 
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'^From the non-linear dynamics point of view, the all zeros 
and all ones are the two stable fixed points of the CA dynam- 
ics, whereas the density equals pc is an unstable fixed point 
of the CA dynamics. 
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in the sequences 11000 and 10011 cannot. Readers can 
readily verify that the total number of cars in the system 
is conserved under the above set of rules. Moreover, one 
can easily convert the above set of traffic rules into a fi- 
nite CA rule table consisting of 2'+^ rules. In addition, 
when pc = 1/2, the modified traffic rule is reduced to 
Wolfram's elementary CA rule 184. 

We define the local car density at each site to be the to- 
tal number of cars contained in that site and in its (7 — 1 
immediate right-neighboring-sites divided by q. Thus, 
the local car density of a site changes when and only 
when a car enters that site, or a car leaves the {q — l)th 
site to its right. For simplicity, a collection of sites is said 
to be in the low density region if the local car density of 
each site in the collection is less than or equal to pc- And 
a high density region is defined as a collection of sites 
with local car densities greater than pc. 

Now, we prove two theorems concerning the distribu- 
tion of high and low density regions under the repeated 
actions of the modified traffic rules. 

Lemma 1 // the initial density p of the system 
is less than or equal to pc, then after at most 
\N{max{q,2p) — l)/<z] + q — 2 timesteps, the local car 
density for every site is less than or equal to pc- 

Proof: Consider a collection of sites with local car density 
greater than pc in the initial system configuration. Using 
the modified traffic rule, in the co-moving frame of a car, 
no site will have car density exceeding {jp -\- 1)/q if that 
car is originally located at a site with local car density 
less than or equal to pc- In addition, such an aggregation 
of car above the density threshold pc will be achieved 
within the first q — 1 timesteps. Since the overall car 
density of the system is less than or equal to pc, so upon 
repeated applications of the modified traffic rules, cars 
will gradually move out of local regions with car density 
exceeding pc (if any). Note that the local car densities 
of those cars moving out of these "high density regions" 
are less than or equal to pc- Moreover, once these cars 
are "dissolved" from a local high density region, the local 
car density of the sites containing these cars will never 
exceed pc unless these cars merge into a high car density 
region in front. 

Suppose none of the cars moving out of a high car den- 
sity region will be stopped by another high car density 
region in front. Then, once all the high car density re- 
gions dissolve completely, no further high density region 
will be formed thereafter. Thus, our assertion that all 
sites will have local car density not exceeding pc asymp- 
totically is true in this case. So, it remains to consider the 
case that some of the cars moving out of a high car den- 
sity region move into (and hence temporarily stopped by) 
a high car density region in front. Since the overall car 
density of the entire system does not exceed pc, at least 
one of the high car density regions will start to dissolve. 
Note that N is finite, and those cars merging into a high 
density region will readily re-dissolve again once they are 



allowed to move. After a finite number of timesteps, all 
the high density regions will disappear. After that, no 
more high density region can form. Thus, our assertion 
that all sites will have local car density not exceeding pc 
is also true is this case. 

Finally, we estimate the number of timesteps required 
to reach this state. Since the updating is taken in parallel, 
the worst case occurs when some of the local high density 
regions are formed during the first q — 2 timesteps. Then 
all the cars merge into a single high car density region 
before they finally dissolve. According to the modified 
traffic rules, if p < q/2, then all cars dissolved from a 
high density region will not be blocked. In this case, it 
takes at most g — 1 timesteps to dissolve p cars from a 
high density region. On the other hand, ii p > g/2, then 
some of the dissolved cars will still be blocked occasion- 
ally. And it takes q — 1 + {2p — q) = 2p — 1 timesteps to 
dissolve p cars from a high density region. Therefore, in 
any case, it requires at most [A^p(max(g, 2p) — l)/p] < 
\N{ma.x{q,2p) — l)/q'] timesteps to completely dissolve 
a high density region making up of iVp cars. Hence, our 
assertion is proved. □ 

Lemma 2 // the initial density p of the sys- 
tem is greater than pc, then after at most 
[A^(g — p)(max(g, 2p) — l)/pg] timesteps, the local car 
density for every site is greater than or equal to pc ■ 

Proof: We may assume that the local car density at each 
site initially is less than or equal to pc- Otherwise, no car 
can move right at the beginning and our assertion is triv- 
ially true. Since the overall car density p is greater than 
Pc all the low car density regions (sites with local car 
density less than pc) are surrounded by high car density 
ones (sites with car density greater than or equal to pc)- 
Thus, cars can gradually move from the high density to 
the low density regions. And using similar arguments as 
in Lemma 1^, after finite number of timesteps, the entire 
system is contained in a single high density region. And 
from that time on, no car in the system can move. 

We move on to estimate the number of timesteps 
required to reach this "frozen" state. Similar to the 
argument in Lemma the worst case occurs when 
there is only one low density region initially. And 
it is easy to verify that the number of timesteps re- 
quired to reach a frozen state for such an initial system 
configuration equals [A^(l — p)(max(g, 2p) — l)/p] < 
\N{q - p){max{q, 2p) - 1)/^] • □ 

Combining Lemmas |l| and H, we conclude that 

Theorem 1 By applying the modified traffic rules for 
[A^(max(g, 2p) — 1) max(q — PTP)lpq] + q — 2 times, an 
initial system configuration will he segregated into one of 
the following three cases: 

(a) If p < Pc, then the local car density at every site is 
less than or equal to pc- In addition, at least one of the 
sites will have local car density strictly less than pc- 
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(b) If p = Pc, then the local car density at every site 
equals pc- 

(c) If p > Pc then the local car density at every site is 
greater than or equal to pc- In addition, at least one of 
the sites will have local car density strictly greater than 

Pc- 

III. MODIFIED MAJORITY RULE 

After segregating the system configuration according 
to its initial density p, the density classification prob- 
lem becomes straight-forward. We consider the follow- 
ing modified majority rule for a given critical density 
Pc = p/q: The state of a site in the next timestep is 
one if there are at least 2p + 1 ones in the 2q + 1 sites 
consisting of itself and the q consecutive left- and right- 
neighboring-sites. Otherwise, the state of this site in the 
next timestep is zero. For example, if pc = 1/2, then the 
states of the middle sites in the next timestep for 10101 
and 01010 are one and zero, respectively. 

Now, we present the results of applying the modified 
majority rule to certain system configurations which are 
of interest. 

Lemma 3 Any system configuration with local car den- 
sity at every site equals pc is a fixed point of the modified 
majority rule dynamics. 

Proof: Consider an arbitrary site a in the system. Since 
the local car density at every site equals pc, there are 
precisely 2p sites in state one among the 2q neighboring 
sites of a. Thus, the majority rule implies that the state 
of site a in the next timestep is equal to its present state. 
Consequently, this system configuration is a fixed point 
of the modified majority rule. □ 

Lemma 4 Suppose the local car density at every site is 
greater than or equal to pc and also that the total density 
of the system is strictly greater than pc- After applying 
the modified traffic rule \N/2{q— 1)] times, the state of 
every site in the system will become one. Similarly, if 
the local car density at every site is less than or equal to 
Pc and also that the total density of the system is strictly 
less than pc . Then after applying the modified traffic rule 
\N/2{q — 1)1 times, the state of every site in the system 
will become zero. 

Proof: For simplicity, we only consider the case that 
p > Pc. The proof for the case that p < pc is similar. 
From Lemma ^ we know that the local car density at 
every site must be greater than or equal to pc after re- 
peated applications of the modified majority rule. More- 
over, if the local car density of a site a exceeds pc, then 
the states of a and its g — 1 left- and right-neighboring- 
sites must be one in the next (and hence all subsequent) 
timestep. In other words, the propagation speed for state 
one is g — 1 sites per timestep both leftward and right- 
ward. Therefore, after \N/2{q — 1)] timesteps, all sites 



in the system will be in state one. Hence, our assertion 
is proved. □ 

Combining Theorem |^, Lemmas ^ and |[ we obtain the 
following CA density classification theorem: 

Theorem 2 (Density Classification By CA) Let 

Pc = p/q be a rational number between zero and one, 
with p, q being relatively prime positive integers. Then 
the density classification problem can be solved using the 
following two CA rules: apply the modified traffic rules 
\N{max{q, 2p) — 1) max(q — p,p)/pq~\ + q — 2 times and 
then followed by the modified majority rule [N/2{q — 1)] 
times. 



IV. DISCUSSIONS 

In summary, we report a way to classify the density 
of ones using two one-dimensional binary CAs provided 
that the density threshold pc is a rational number. Our 
result, therefore, generalizes that of Fuks 0. Besides, 
our construction also takes care of the case when the 
density p is equal to its critical value pc — something 
Fuks does not consider seriously. Because of the local 
nature of CA rules, we believe that the time complex- 
ity of a general density classification problem using any 
combination of CA rules is at least 0{N). Thus, apart 
from a constant speedup, our set of density classification 
CA rules is probably the least time-consuming. 

As we have discussed in Sec. |, if the number of site 
in the system N is finite, it suffices to restrict ourselves 
to consider rational density classification. Any irrational 
density can be well approximated by a corresponding ra- 
tional density. Actually, to approximate an irrational 
density by some better and better rational numbers, 
the size of the CA rule table grows. Thus, we believe 
that combinations of CA rules is unlikely to be powerful 
enough to classify irrational density in an infinite one- 
dimensional system. 
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